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ABSTRACT 



In order to understand more fully the wind generation of 
water waves, the oceanographer should have a better knowledge 
of the turbulence structure of the wind field immediately above 
a water surface. 

Measurements of atmospheric turbulence about three inches 
above the surface of a lake, using a hot-wire anemometer, were 
made. The analog signals were recorded as DC voltages on a 
magnetic tape recorder. The analog data were converted to 
digital form to permit subsequent computation of power spectral 
densities using a high-speed digital computer. 

Average power spectral densities from 0 to 20, 0 to 500, 
and 0 to 1000 cps, respectively, are presented. 
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1 • Introduction 



The fact that the wind generates and maintains waves on the 
ocean surface and provides the driving force for many ocean 
currents has been well known for many years. The exact mechanism 
by which the wind produces these effects is, however, not well 
known. Prior to understanding completely the nature of the air- 
sea interaction, the oceanographer must describe accurately the 
properties of the wind field immediately over the ocean surface. 

Most analytical calculations of the rate of energy transfer 
from the wind to progressive water waves have used the assumption 
that the air flow over the water surface is basically laminar. 

This assumption, which, in effect, neglects atmospheric turbu- 
lence, is not realistic if the wave-generation process is to be 
completely understood. 

Measurements of turbulence in the wind field in the first 
few inches above a water surface would contribute to a better 
understanding of the interactions between air and sea. The object- 
ives of this research, therefore, have been threefold: 

1. To develop and test a data-collection system capable of 
measuring and recording turbulent fluctuations in the natural 
wind immediately over a water surface. 

2. To make recordings of wind turbulence in the first few 
inches above a water surface. 

3. To compute power spectral densities (PSD*s) from the 
recorded data. 

Although the system developed is restricted to use in a relatively 
shallow fresh-water lake, the results obtained and the experience 
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gained may serve as a foundation for further work in this 
field at the Naval Postgraduate School (NPGS), 

2 . Description of Observations 

In order to measure and record the turbulent fluctuations 
in the downstream component of the wind field, it was necessary 
to devise a data-collection system suited to the task. Following 
a careful search of the literature pertaining to this problem, 
and after a survey of equipment available at the NPGS, a basic data- 
collection system was selected. Figure 1 is a block diagram of the 
system used in this research. 

The instrument used for measuring the turbulence in the wind 
field was a constant temperature hot-wire anemometer. Security 
Associates Model 200. This anemometer, which consists of a high 
gain constant temperature amplifier coupled to an analog computer 
system, gives an output which is linear in velocity and zero for 
zero velocity. The accuracy of the analog transfer function is 
better than 0.5 percent over the operating range (DC-50 kilocycles 
per second). The hot-wire probes consisted of tungsten wires about 
0.0015 inches in diameter and 0.090 inches in length. 

The signal from the anemometer was recorded on an FM channel 
of a magnetic tape recorder, the AMPEX Model CP- 100: a 14 channel, 

solid state, portable instrument. Recordings were made at a tape 
speed of 3.75 inches per second, giving a frequency bandwidth from 
0 to 1250 cycles per second (cps). Before recording, the recorder 
was calibrated in accordance with the standard procedure outlined 
in the manufacturer’s instruction book. The signal from the 
anemometer was biased so that the turbulent fluctuations were 
recorded about zero volts DC. This was necessary to insure that 
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the recorder signal remained within the voltage limits of the 
analog-to-digital (A/D) converter used in data processing. The 
biased signal was then amplified (x3) and monitored on an oscil- 
loscope to insure that the correct biasing was maintained. 

Before collecting data in the field, the system was tested 
using fan-generated wind over a wave tank. The test signals 
were recorded and analyzed. The results of the test runs indicated 
that the system was functioning properly. 

In order to maintain the hot-wire probe at a constant distance 
above the mean water level and normal to the mean wind field, a 
special wind vane and probe mount were constructed (Figure 2). 

At 1600, 7 April 1967, the equipment was set up at Roberts 
Lake, Seaside, California. This site was selected because of its 
relatively unbroken wind field, its close proximity to the 
Naval Postgraduate School, and the availability of electric power 
for the equipment. Roberts Lake lies 500 yards inland from the 
Seaside Beach, which borders Monterey Bay. The lake, which is 
about 500 yards long and 200 yards wide at the center, is approx- 
imately oval shaped, with the long axis running parallel to the 
beach. 

The wind vane, with the hot-wire probe mounted about 3 inches 
above the mean water level, was set in about 2.5 feet of water at 
the downwind end of the lake. This position gave an unbroken 
fetch of approximately 75 yards. The signal was recorded for 
30 minutes with the wind vane free to fluctuate in the wind field. 
During this period, it was noted that the vane had a tendency to 
oscillate about the direction of the mean wind. Upon completion 
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of the 30 minute recording period, the vane was fixed so that 
the hot-wire probe was approximately normal to the mean wind 
direction; the signal was then recorded for an additional 10 
minutes . 

During the recording period, the average wind speed varied 
between 6 and 8 knots, with the gusts up to 12 knots. The signif- 
icant wave height was 1 to 2 inches, the significant wave length 
15 to 20 inches, and the significant wave period 1 to 2 seconds. 
Since no accurate measuring devices were available, these values 
were subjectively estimated. 

3. Method of Analysis 

Based upon the results of Pond et^ ,al. (1966), it was expected 
that a large percentage of the turbulent kinetic energy would lie 
in the frequency range from 0 to 5 cps. The analog equipment 
available at the NPGS was not capable of meaningful analysis in this 
frequency range. For this reason a digital method of analysis was 
chosen. 

Atmospheric turbulence is believed to approximate a stationary, 
random process. It is suitable, therefore, to apply the Blackman 
and Tukey (1958) method of analysis to compute the PSD. This method 
uses the fact that the PSD can be expressed as the Fourier transform 
of the autocorrelation function. The method readily lends itself 
to the computation of the PSD using a high-speed digital computer. 

In order to convert the continuous record of analog data, 
which were in the form of varying DC voltage levels on a magnetic 
tape, to one of discrete data suitable for digital computation, it 
was necessary to sample the data at equally spaced time intervals. 
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This was accomplished by using a CDC 160 Computer coupled to an 
A/D converter for input and CDC 163 Tape Unit for output. The 
digitizing program, PROGRAM POWSPEC, is capable of sampling the 
data at a maximum rate of 5000 discrete sample points per second. 

This allows computation of the PSD*s in the frequency range from 
0 to 2500 cps. 

The final step in the analysis of the data was the computation 
of the PSD's, using PROGRAM POWSPEC, a FORTRAN 63 program, in a 
CDC 1604 Computer. Detailed discussions of the A/D conversion 
and PSD program are given in Appendices A and B. 

4. Results and Conclusions 

The record was processed to obtain estimates of the power 
spectral densities from 0 to 1000 cps, 0 to 500 cps, and 0 to 20 
cps, using low-pass filters to filter out undesirable high-frequency 
noise and to prevent "aliasing" of high-frequency energy into the 
low-frequency end of the spectrum. The record was processed also 
from 0 to 1000 cps and 0 to 500 cps without filters for comparison 
with the filtered data. The power spectra were further categorized 
according to whether the wind vane was free to oscillate about the 
mean wind or was held fixed into the mean wind. Due to the unusual 
design of the hot-wire probe mount, calibration of the anemometer 
system was impractical. Only the relative contributions of spectral 
energy from the various frequencies, therefore, were measured. 

Five spectral estimates in each category were averaged to obtain 
an approximation of the power spectrum in that category. At the 80- 
percent confidence level, an estimate of any single value in the 
average spectra lies between the percentages of the actual value 
as tabulated below: 
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Frequency 
Range (cps) 


Number of 
Data Points 


Number of 
Lags 


Percent of 
Actual Value 


0-20 


4000 


500 


79-123 


0-500 


4000 


125 


89-111 


0-1000 


4000 


250 


85-115 . 



The average power spectra are graphically presented in Figures 3 
through 9. The values from which these spectra were plotted are 
included as Tables IV through X in Appendix C. 

An examination of the 0 to 1000 cps and 0 to 500 cps spectra 
(Figures 5 through 9) shows that most of the energy lies in the 
frequency range from 0 to 50 cps. Energy peaks occurring beyond 
about 50 cps appear to be randomly distributed with frequency. 

Some possible causes of these peaks are mechanical vibrations in 
the wind vane, noise in the electronic equipment in the data-col- 
lection and processing systems, and changes in the level of turbu- 
lence in the wind field. The addition of filters during data 
processing produced no meaningful changes in the average spectra. 

A high-resolution analysis of the 0 to 20 cps frequency band, 
using a 20 cps low-pass filter, was made. The average spectra in 
this frequency band, with the vane free to oscillate and with the 
vane fixed, are shown in Figures 3 and 4, respectively. Although 
the total energy computed in the fixed-vane case is greater, the 
differences between the two spectra are not considered significant. 

The significant frequency range of the water waves present 
during the observations, based on a subjective estimate of the 
significant wave lengths, was computed using the small-amplitude 
relation of frequency to wave length (f = [5.12/L]^ 2 ). The 
frequency range of the water waves thus computed was 1.7 to 2.1 cps 
this range has been indicated by the vertical dashed lines on 
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Figures 3 and 4. On each of these two average spectra, an energy 
peak appears within the frequency range of the water waves. These 
peaks are attributed to the influence of the water-wave undulance 
on the turbulence structure of the wind field. It should be noted 
that, although the probe was only about twice as high above the 
mean water surface as were the water waves, the undulance of the 
water waves seems to contribute a relatively small percentage of 
the total turbulent energy. 

PROGRAM POWSPEC estimates the total energy in two different 
ways: (1) by computing the variance of the 4000 data points and 

(2) by integrating the power spectrum by application of the 
trapezoidal rule. The ratio of the two energy values thus obtain- 
ed gives the fraction of the total energy accounted for in the 
power spectrum. The value of this ratio is called XFACT in the 
program. If all of the energy is accounted for in the power 
spectrum, XFACT will equal unity. 

The average XFACT in the 0 to 500 cps and the 0 to 1000 cps 
power spectra, using a 4 cps frequency interval, was about 0.80, 
indicating that the power spectral estimates did not account for 
nearly 20 percent of the energy. Spectra in these frequency ranges 
were recomputed using a 1 cps frequency interval, and XFACT increased 
to about 0.90. The shapes of the spectra were unchanged. The power 
spectral estimates for the 0 to 20 cps spectra, using a 0.1 cps 
frequency interval, accounted for about 94 percent of the energy. 
XFACT apparently can be made to approach unity by decreasing the 
frequency interval. This is probably true because as the number of 
estimates is increased the trapezoidal rule provides a more accurate 
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approximation of the value of the integral of the power spectrum. 

The confidence in any single value in the spectrum, however, is 
decreased when the number of estimates is increased. 

In every power spectrum computed a large spectral density 
value was obtained at zero frequency. In every average spectrum, 
the maximum spectral density value was obtained at one frequency 
interval abo^ve zero frequency. Precautions are taken in PROGRAM 
POWSPEC to prevent unrealistic values of power density at or near 
zero frequency in that the mean of the data is subtracted from 
each data point before computing the PSD's, as suggested by 
Blackman and Tukey (1958). The significance of the large spectral 
density values remaining at the very low end of each spectrum is 
not f u 1 ly under s t ood . 

5 ♦ Recommendations for Future Work 

1. A detailed analysis of the very low frequency end of the 
spectrum should be made to determine more accurately the distribution 
of energy in this frequency range. 

2. A float-mounted probe should be used to measure the turbu- 
lence in the wind field at a constant distance above the actual 
water surface. A comparison between spectra thus obtained and the 
spectra obtained at a fixed level above the mean water surface should 
be made. 

3. Crossed hot-wire arrays anda multichannel hot-wire anemometer 
should be used to measure the downstream and transverse components 

of atmospheric turbulence in order to compute co-spectral densities 
and Reynolds stresses. 
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FIGURE-1 



BLOCK DIAGRAM OF DATA COLLECTION SYSTEM 
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FIGURE 2. WliJD VANE Al.D PROBE MOUNT 
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FIGURE 3. AVERAGE POWER SPECTRUM, 0-20 ops, FILTER, VANE FREE 
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FIGURE 5* AVERAGE POWER SPECTRUM, 0-500 ops, NO FILTER VANE FREE 
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Frequency ( ops ) 

FIGURE 6. AVERAGE POWER SPECTRUM, 0-500 cps, NO FILTER, VANE FIXED 
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PIGURE «. AVERAGE POWER SPECTRUM, 0-500 cps, FILTER, VANE PIXED 
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FIGURE 9 . AVERAGE POWER SPECTRUM, 0-1000 ops, FILTER, VANE PIXED 
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APPENDIX A 



ANALOG TO DIGITAL DATA CONVERSION 
1 . Introduction 

The analog data, recorded on magnetic tape as a varying 
DC voltage, is converted to digital form in order to use a high- 
speed digital computer for analysis. The output data from the 
AMPEX CP- 100 Recorder is fed through an analog- to-digital (A/D) 
converter. The digitizing computer, CDC Model 160, alternately 
activates the Ayb converter for data input and the tape unit, 

CDC Model 163, for data output. 

The CDC 160 Computer has a 4096-cell memory. Each cell 
contains 12 bits. The top 96 cells are used for programming, 
while the lower 4000 cells are reserved for data storage. A 
block diagram of the digitizing system is shown in Figure A-l. 

When the A/D converter is activated and the input instruction 
executed, the analog input data from the recorder is sampled and 
a 12-bit number, representing signal amplitude, is sent to the 
CDC 160 and stored. This process is repeated in rapid succession. 
Blocks of up to 4000 samples can be stored and the maximum sampling 
rate is 5000 samples per second. 

The normal mode of operation is to initiate each block in a 
run by a manual gating pulse on the A/D converter. The computer 
controls the inter-sample delay time. Automatic pulsing, in 
which the computer controls the gating of each block in a run, is 
also available. 
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2. Calibration 



Prior to the actual digitizing process, the system must 
be calibrated using a test program, TEST 160 (Table I), In 
PROGRAM TEST 160, cells 0000 to 0017 contain a program which 
accurately sets the -5.0 volt bias voltage required by the 
A/D converter. If the bias voltage is properly set at -5.0 
volts, with no input signal to the A/D converter, the value 
of the conversion, as indicated in cell 0004, should be 0000 
or 7777. Small variations in the bias voltage, which are to 
be expected during the course of a run, are eliminated by 
removing the mean of the data during computation of the 
power spectral density. 

Cells 0020 to 0045 in PROGRAM TEST 160 contain a program 
to set the inter-sample delay time accurately. The input 
request pulses to the A/D converter are observed on an oscil- 
loscope. The inter-sample delay time can be adjusted to the 
desired value by manually changing the timer-control word in 
the program (cell 0042). Figure A-2 is a plot of timer-control 
wbrd (octal) versus inter-sample delay time. 

3 . Digitizing the Analog Data 

After the calibration procedure is completed, PROGRAM 



DIGITIZE (Table 


II) 


is read into the computer. The following 


manual entries 


into 


PROGRAM DIGITIZE are required 


: 


Cell No. 




Contents 


Remarks 


0066 




Timer Control Word 


See Figure A-2 


0067 




Spare Identifier 


Any number ^7777 


0070 




Initial Run Number 


Set as 0001 


0071 




Number of runs desired 




0072 




Number of samples/block 


^7640 g = 4000 10 


0073 




Number of blocks/run. 






Each block of digitized data is labeled with a 16-digit 
octal identifying number so that the block can be recalled 
for computation in the main PSD program. Four quantities are 
specified in the identifier: 

1. Current run number - steps once per run 

2. Spare ID - may be set to any number ^7777g 

3. Current block number in the run 

4. Number of blocks per run. 

As an example, assume that ten runs, forty blocks per run, 
with the spare ID set at 0015, are being used. If the fourth 
and fifth blocks of the sixth run are recalled, the identifiers 
would appear in the CDC 1604 as follows: 

0006 0015 0004 0050 

0006 0015 0005 0050. 

Exact records of the identifiers must be kept while digitizing, 
as this is the only means of recalling the data desired for 
computation in the main program. 
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FIGURE A-l 



ANALOG TO DIGITAL CONVERSION SI3TEM 
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A -2 . INTER-SAMPLE DELAY TIME VS TIMER CONTROL WORD 
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TABLE I 



PROGRAM TEST 160 



Cell 


Contents 


Code 


0000 


7500 


EXFOO 


1 


2410 


2410 


2 


7500 


EXFOO 


3 


1401 


1401 


4 


7600 


1NA 


5 


4070 


STD 70 


6 


7500 


EXFOO 


7 


2401 


2401 


0010 


7303 


OUT 03 


11 


0072 


0072 


12 


6102 


NZF02 


13 


0070 


0070 


14 


0400 


LDN00 


15 


4071 


STD71 


16 


6414 


ZJB14 


17 


7700 


HLT00 


0020 


7500 


EXFOO 


21 


2410 


2410 


22 


7500 


EXFOO 


23 


1401 


1401 


24 


7600 


INA 


25 


4143 


STI43 


26 


0300 


NOP 


27 


0300 


NOP 


0030 


0300 


NOP 


31 


0300 


NOP 


32 


0300 


NOP 


33 


0300 


NOP 


34 


0300 


NOP 


35 


0300 


NOP 


36 


2042 


LDD42 


37 


0701 


SBN01 


0040 


6501 


NXB01 


41 


7044 


JPI44 


42 


M 


M 


43 


0043 


0043 


44 


0022 


0022 


45 


7701 


HLT 01 



Explanations 

Begin Bias-set 
Set Enable 

Call A/D Channel 1 
Input to "A" in Reg. 
Store "A” in 0070 

Call D/A Channel 1 

Output Cell 0070 to 
Cell 0071 INC. 



Load zero's. 

Store in 0071 
Jump back to 0000 
HLT. End Bias-set 

Begin Timer 
Set Enable 
Call A/D Ch. 1 



Store Sample in 0043 
No Operation-NOP 
These are Time 

Dummies to 

Match Program DIGITIZE 
Timing 



Load Timer Control Word 
Subtract i 

If not zero, go back 1 
If zero, lump to 0022 
Timfcr Control Word 



End Timer 



Note: M is manual entry. 
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TABLE II 



PROGRAM DIGITIZE 



Cell 


Contents 


Code 


0000 


0101 


PTA 


1 


0603 


ADN03 


2 


7064 


JPI64 


3 


7500 


EXF00 


4 


2410 


2410 


5 


7500 


EXF00 


6 


1401 


1401 


7 


7600 


INA 


0010 


4176 


STI76 


11 


2076 


LDD76 


12 


3465 


SBN65 


13 


6134 


NZF34 


14 


7500 


EXF00 


15 


2400 


2400 


16 


2074 


LDD74 


17 


4160 


STI60 


0020 


2067 


LDD67 


21 


4161 


STI61 


22 


2075 


LCD 7 5 


23 


4162 


STI62 


24 


2073 


LDD73 


25 


4163 


STI63 


26 


7500 


EXF00 


27 


2111 


2111 


0030 


7303 


OUT 03 


31 


0000 


C 


32 


6102 


NZF02 


33 


0133 


0133 


34 


2246 


LDF46 


35 


4076 


STD76 


36 


2075 


LDD75 


37 


3473 


SBD73 


0040 


6155 


NZF55 


41 


0401 


LDN01 


42 


4075 


STD75 


43 


5455 


A0D55 


44 


3454 


SBD54 


45 


6153 


NZF53 


46 


6061 


ZJF61 


47 


5476 


AOD76 



Explanations 



Jump to INITIAL 

Set Enable 

Call A/D Channel 1 
INPUT 

Store sample in (0076) 

Enough samples yet? 

If not, go to cell 0047 

Clear Enable 

Load Current Run No. 

Store in cell 0133 

Load Spare ID. 

Store in cell 0134 
Load Current Block No. 

Store in cell 0135 
Load No. of Blocks/Run 
Store in cell 0136 

Call 163 M.T.U. 

Output from 
computed L.W.A. 
to 

0133 

Set A-0137 

Reset running storage address 
Enough Blocks yet? 

If not, go to cell 0115 

Reset block no. 

Update total block count 
Check if capacity exceeded 
If not, go to cell 0120 
If not, go to cell 0127 
Begin intersample delay loop 
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Cell 

0050 

51 

52 

53 

54 

55 

56 

57 

0060 

61 

62 

63 

64 

65 

66 

67 

0070 

71 

72 

73 

74 

75 

76 

77 

0100 

101 

102 

103 

104 

105 

106 

107 

0110 

11 

12 

13 

14 

15 

16 

17 



TABLE xx (continued) 



Contents Code 

2066 LDD66 

0701 SBN01 

6501 NZB01 

7056 JPI56 

0760 0760 

0000 C 

0005 0005 

0003 0003 



0133 0133 

0134 0134 

0135 0135 

0136 0136 

0100 0100 

C C 



M M 

M M 



M 

M 

M 

M 

C 

C 

C 

C 



M 

M 

M 

M 

C 

C 

c 

c 



'1-077 STD77 

2200 LDCOO 

0137 0137 

4076 STD 7 6 

3072 ADD72 

4031 STD 31 

0701 SBN01 

4065 STD65 



2070 LDD70 

4074 STD 7 4 

0401 LDN01 

4075 STD75 

7077 JPI77 

5475 A0D75 

5455 AOD55 

7057 JPI57 



Explanations 

Load intersample delay word 
If not zero, go back 1 
If zero, go to cell 0005 
0760-M.T.U. block capacity 
Total No. Blocks written 
Address 
Address 



Address of Current Run No. ID 
Address of Spare ID of ID 
Address of Cur. Block No. ID 
Address of Tot. Blocks ID 
Address of INITIAL 
Address of last word of data 
Intersample delay word 
Spare ID anything 

Initial Run No. Set 1 
No. runs desired 
No. samples/block 7640g 
No. blocks /run 
Current run no. 

Current Block Number 
Running Storage Address 
JUMP CONTROL 

BEGIN INITIAL 
Set "A" -0137 

Initialize running storage address 
Compute last address 
Store in 0031 

Store in 0065 



Initialize Run No. 

Initialize Block No. 

END INITIAL 

Update Current Block No. 
Update Current Block Count 
Go to 0003 
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TABLE II (continued) 



Cell 


Contents 


Code 


0120 


2074 


LDD74 


21 


3471 


SBD71 


22 


6103 


NZF03 


23 


2074 


LDD74 


24 


7701 


HLT01 


25 


5474 


AOD74 


26 


7057 


JPI57 


27 


7500 


EXF00 


0130 


1111 


1111 


31 


2055 


LDD55 


32 


7702 


HLT02 


33 


C 


C 


34 


C 


C 


35 


C 


C 


36 


C 


C 


0137 


D 


D 




D 


D 


7776 


D 


D 



Have enough runs 
been done yet? 

If no, go to 0125 

If yes, display last run no. 

Halt 01 

Update current run no. 

Go to 0003 
Call E.O.F. 

Write an EOF Mark 
Load block count 
Halt 02 
RUN #IB 
SPARE #ID 
BLOCK #I.D 
BLOCKS /RUN ID. 

Data Storage 

Data Storage 



Note : 

C is PROGRAM COMPUTED ENTRY 
M is MANUAL ENTRY 
D is DATA STORAGE 
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APPENDIX B 



PROGRAM FOR COMPUTING POWER SPECTRAL DENSITIES 

1 . Introduction 

After the digitizing is completed, the magnetic tape 
contains binary numbers, the form of which is suitable for 
the computation of power spectral density (PSD) by a high- 
speed digital computer. PROGRAM POWSPEC, a FORTRAN 63 program 
which uses the Tukey method with a banning lag window for 
smoothing, is used to compute the PSD (Table III). This 
program is available at the Naval Postgraduate School. 

2. Data Recall 

The data in digital form are recalled in PROGRAM POWSPEC 
by calling SUBROUTINE DATA, which uses two subroutines, FINDIT 
and UNPACK. SUBROUTINE FINDIT locates the data on the tape and 
SUBROUTINE UNPA.C I unpacks the block of numbers into the memory 
of the CDC 1604. The data to be recalled is designated by its 
16-digit octal identifying number. 

3 . Main Program 

The unpacked data is in the memory of the CDC 1604 in 
integer format. PROGRAM POWSPEC converts the data to floating 
point format, divides each data point by -409.6 (which reconverts 
the data point to the voltage which it represents), and removes 
the mean of the data. The program then computes the mean square 
and the autocorrelation for each lag. By means of a Fourier 
cosine transform the power spectrum is computed from the auto- 
correlations. The PSD with respect to frequency is then normalized 
with respect to the energy at zero frequency. The program prints 
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out the autocorrelations and the PSD's versus the frequencies, 
and plots a normalized frequency spectrum. 

4. Data Cards 

In order to compute the frequency spectrum certain data 
cards are required by the main program. These cards are as 
follows : 

1. Data Control Card 

a) ITAPE - indicates whether data is read in from 
tape or cards. ITAPE = 0 for tape data and 1 
for card data. 

b) IDELTA - indicates that every IDELTA data point 
is to be used in the computation of the frequency 
spectrum. If every data point is to be used, 

IDELTA « 1. 

c) DELTAT - the inter-sample delay time used when 
digitizing data. 

d) ISTART - indicates which data point is the starting 
point for the computation. If the computations 
begin with the first data point, ISTART = 1. 

e) IMAX - the maximum number of data points to be used 
for computation. If data is read in from tape, 
using the CDC 160 digitizing process, IMAX ^4000. 

f) M - the number of estimates to be used in computing 
the frequency spectrum. 



37 . 



2. Graph Control Card 



a) ABS - the number of units per inch for the 
abscissa. This number is rounded to one signif- 
icant figure by the program. 

b) ORD - the number of units per inch for the 
ordinate. This number is rounded to one signif- 
icant figure by the program. 

c) IWD - width of graph in inches (1 ^ IWD ^ 9). 

d) IHI - height of the graph in inches (1 ^ IHI ^ 15). 

3. Other Data Cards 

a) NUMHDR - the number of blocks to be computed. 

b) Block Identifier - a 16-digit octal number 
assigned during the digitizing process. 

c) Graph Identifier - an eight character alpha- 
numeric graph identifier. 

5 . General Comments 

The maximum frequency (f max ) to which the frequency spectrum 
may be computed is controlled by the selection of DELTAT (AT). 

This must be determined before digitizing. The following equation 
applies : 



f 

max 



i / [ 2 at ] . 



The number of intervals between zero frequency and the 
maximum frequency is determined by M, the number of estimates, 
in accordance with the following equation: 

Af = 1/[2ATM]. 
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GO TO 30 

START TAPE DATA READ IN— USE THIS BRANCH FOR PRODUCTION 
NUMHDR=NO. OF HEADINGS 
READ 1 10 » NUMHDR 
DO 7 N*l. NUMHDR 
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TABLE III (continued) 
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TABLE III (continued) 
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APPENDIX C 



TABLES OF COMPUTED AVERAGE SPECTRAL DENSITIES 
Five power spectra were averaged for each average power 
spectrum shown in Figures 3 through 9. The average values of 
power spectral density versus frequency from which the average 
spectra were drawn are listed in Tables IV through X. 
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TABLE IV 





AVERAGE POWER 


SPECTRAL DENSITIES 






0-20 cps, FILTER, VANE FREE 




Frequency 


(cps) Energy 


Frequency (cps) 


Energy 


0.0 


1.02972 


4.5 


.00628 


0.1 


2.40273 


4.6 


.00524 


0.2 


2.17715 


4.7 


.00656 


0.3 


1.56216 


4.8 


.00739 


0.4 


.86354 


4.9 


.00734 


0.5 


.41093 


5.0 


.00758 


0.6 


.24272 


5.2 


.00507 


0.7 


.19241 


5.4 


.00462 


0.8 


.14300 


5.6 


.00503 


0.9 


.12687 


5.8 


.00565 


1.0 


.09937 


6.0 


.00358 


1.1 


.07135 


6.2 


.00352 


1.2 


.05696 


6.4 


.00255 


1.3 


.04846 


6.6 


.00329 


1.4 


. 04885 


6.8 


.00366 


1.5 


.04625 


7.0 


.00357 


1.6 


.03681 


7.2 


.00288 


1.7 


.03397 


7.4 


.00228 


1.8 


.03451 


7.6 


.00248 


1.9 


.03137 


7.8 


.00312 


2.0 


.03408 


8.0 


.00238 


2.1 


.03898 


8.2 


.00247 


2.2 


.03691 


8.4 


.00175 


2.3 


.03157 


8.6 


.00161 


2.4 


.02312 


8.8 


.00130 


2.5 


.01920 


9.0 


.00117 


2.6 


.02404 


9.2 


.00126 


2.7 


.02467 


9.4 


.00096 


2.8 


.01899 


9.6 


.00111 


2.9 


.01789 


9.8 


.00100 


3.0 


.01956 


10.0 


.00070 


3.1 


.01738 


11.0 


. 00084 


3.2 


.01548 


12.0 


.00038 


3.3 


.01781 


13.0 


.00045 


3.4 


.01733 


14.0 


.00040 


3.5 


.01250 


15.0 


.00024 


3.6 


.01008 


16.0 


.00017 


3.7 


.01171 


17.0 


.00011 


3.8 


.01288 


18.0 


.00010 


3.9 


.01162 


20.0 


.00006 


4.0 


.01053 






4.1 


.00980 


Average XFACT = 


.95139 


4.2 


.00891 


Average A(0) = 


.25695 


4.3 


.00885 






4.4 


.00849 
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TABLE V 



AVERAGE POWER SPECTRAL DENSITIES 
0-20 cps, FILTER, VANE FIXED 



Frequency (cps) 


Energy 


?' Frequency(cps) 


Energy 


0.0 


1.58421 


4.5 


.00827 


0.1 


2.74913 


4.6 


.00884 


0.2 


1.87637 


4.7 


.00917 


0.3 


1.06050 


4.8 


.00836 


0.4 


.56193 


4.9 


.00744 


0.5 


.39477 


5.0 


.00696 


0.6 


.34068 


5.2 


.00834 


0.7 


,26455 


5.4 


.00651 


0.8 


.21270 


5.6 


.00665 


0.9 


.15058 


5.8 


. 0047 1 


1.0 


.10860 


6.0 


.00497 


1.1 


.10138 


6.2 


.00661 


1.2 


.08719 


6.4 


.00438 


1.3 


.08060 


6.6 


.00516 


1.4 


.08737 


6.8 


. 00494 


1.5 


.07508 


7.0 


.00383 


1.6 


.05162 


7.2 


.00534 


1.7 


.04714 


7.4 


.00431 


1.8 


.04651 


7.6 


.00357 


1.9 


.04443 


7.8 


.00452 


2.0 


.04520 


8.0 


.00395 


2.1 


.03793 


8.2 


.00333 


2.2 


.03017 


8.4 


.00285 


2.3 


.03243 


8.6 


.00314 


2.4 


.03381 


8.8 


.00209 


2.5 


.02888 


9.0 


.00221 


2.6 


.02345 


9.2 


.00174 


2.7 


.02105 


9.4 


.00186 


2.8 


.02309 


9.6 


.00195 


2.9 


.02518 


9.8 


.00166 


3.0 


.02680 


10.0 


.00135 


3.1 


.02716 


11.0 


.00141 


3.2 


.02415 


12.0 


.00091 


3.3 


.02163 


13.0 


.00058 


3.4 


.01860 


14.0 


. 00038 


3.5 


.01473 


15.0 


.00036 


3.6 


.01297 


16.0 


. 00022 


3.7 


.01286 


17.0 


.00018 


3.8 


.01403 


18.0 


.00019 


3.9 


.01516 


20.0 


.00007 


4.0 


.01443 






4.1 


.01306 


Average XFACT = 


.92798 


4.2 


.01136 


Average A(0) = 


.27124 


4.3 


.00940 






4.4 


.00853 
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TABLE VI 



AVERAGE POWER SPECTRAL DENSITIES 
0-500 cps, NO FILTER, VANE FREE 



Frequency (cps) 


Energy 


Frequency (cps) 


Energy 


0 


.12076 


188 


.00056 


4 


.12954 


192 


00027 


8 


.01088 


212 


.00034 


12 


.00348 


232 


.00041 


16 


.00227 


252 


.00024 


20 


.00159 


272 


.00017 


24 


.00120 


292 


00020 


28 


.00098 


312 


.00020 


32 


.00085 


332 


.00038 


36 


.00069 


352 


.00017 


40 


.00054 


372 


.00019 


44 


.00045 


392 


00047 


48 


.00045 


412 


.00028 


52 


.00047 


420 


00038 


72 


.00034 


424 


.00050 


92 


.00045 


428 


.00053 


104 


.00057 


432 


.00046 


112 


.00040 


452 


.00026 


132 


.00047 


472 


00042 


152 


.00021 


492 


00031 


172 


.00018 


496 


. 00040 


180 


.00048 


500 


.00023 


184 


.00084 







Average XFACT = .80523 
Average A(0) = .24117 
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TABLE VII 



AVERAGE POWER SPECTRAL DENSITIES 
0-500 cps, NO FILTER, VANE FIXED 



Frequency (cps) 


Energy 


Frequency (cps) 


Energy 


0 


.11824 


140 


.00065 


4 


.12909 


152 


.00017 


8 


.01402 


172 


.00019 


12 


.00562 


192 


.00013 


16 


.00394 


212 


.00012 


20 


.00261 


232 


.00025 


24 


.00193 


252 


.00014 


28 


.00159 


272 


.00030 


32 


.00145 


292 


.00010 


36 


.00121 


312 


. 00004 


40 


.00098 


332 


.00009 


44 


.00086 


352 


00005 


48 


.00076 


372 


.00003 


52 


.00068 


392 


.00009 


64 


.00083 


404 


.00022 


68 


.00086 


408 


.00014 


72 


.00058 


412 


.00007 


92 


.00061 


432 


.00007 


100 


.00044 


448 


.00057 


104 


.00039 


452 


.00057 


112 


.00027 


472 


.00055 


132 


.00113 


492 


.00007 


136 


.00129 


500 


.00002 



Average XFACT = .80876 
Average A(0) = .28710 
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TABLE VIII 



AVERAGE POWER SPECTRAL DENSITIES 
0-500 cps, FILTER, VANE FREE 



Frequency (cps) 


Energy 


Frequency (cps) 


Energy 


0 


.13247 


212 


.00021 


4 


.14017 


232 


.00028 


8 


.01024 


236 


.00064 


12 


.00390 


240 


.00084 


16 


.00220 


244 


.00048 


20 


.00154 


248 


.00013 


24 


.00118 


252 


.00009 


28 


.00088 


272 


.00005 


32 


.00073 


292 


.00003 


36 


.00063 


312 


.00005 


40 


.00055 


320 


.00022 


44 


.00045 


332 


.00006 


48 


.00043 


352 


.00008 


52 


.00041 


372 


.00006 


72 


.00025 


392 


.00012 


80 


.00049 


412 


.00013 


84 


.00031 


432 


.00021 


92 


.00018 


436 


.00029 


112 


.00018 


440 


.00026 


132 


.00030 


452 


.00013 


152 


.00012 


472 


.00020 


172 


.00015 


492 


.00016 


192 


.00017 


500 


.00008 



Average XFACT = .79104 
Average A(0) = .25916 
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TABLE IX 



AVERAGE POWER SPECTRAL DENSITIES 
0-500 cps, FILTER, VANE FIXED 



Frequency (cps) Energy 



0 


. 10440 


4 


.11572 


8 


.01525 


12 


.00688 


16 


.00474 


20 


.00296 


24 


.00210 


28 


.00190 


32 


.00172 


36 


.00136 


40 


.00116 


44 


.00119 


48 


.00123 


52 


.00106 


72 


.00063 


92 


.00043 


112 


.00059 


132 


.00071 


152 


.00028 


172 


.00029 


192 


.00030 


204 


.00042 


208 


.00038 


212 


.00039 



Frequency (cps) Energy 



232 .00036 
252 .00027 
272 00037 
292 .00029 
312 .00020 
316 .00022 
332 .00030 
352 .00026 
372 .00023 
392 .00040 
396 .00043 
400 .00041 
412 .00029 
432 .00018 
436 .00035 
440 .00065 
444 .00066 
452 .00040 
472 .00075 
476 .00088 
480 .00046 
492 .00012 
500 .00007 



Average XFACT = . 82849 
Average A(0) = .27232 
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TABLE X 



AVERAGE POWER SPECTRAL DENSITIES 
0-1000 cps, FILTER, VANE FIXED 



Frequency (cps) 


Energy 


Frequency (cps) 


Energy 


0 


.09235 


292 


.00005 


4 


.11094 


312 


.00007 


8 


.02349 


332 


.00006 


12 


.00879 


352 


.00007 


16 


.00638 


372 


.00005 


20 


.00488 


392 


.00006 


24 


.00298 


412 


.00011 


28 


.00203 


432 


.00003 


32 


.00166 


452 


.00006 


36 


.00159 


472 


.00005 


40 


.00146 


492 


.00007 


44 


.00135 


512 


.00027 


48 


.00120 


552 


.00001 


52 


.00111 


592 


.00001 


72 


.00087 


632 


.00002 


92 


.00048 


672 


.00003 


112 


.00068 


712 


.00003 


132 


.00053 


752 


.00010 


152 


.00035 


792 


.00016 


172 


.00022 


832 


.00010 


192 


.00014 


872 


.00017 


212 


.00013 


912 


.00017 


232 


.00028 


952 


.00008 


252 


.00008 


992 


.00012 


272 


.00007 


1000 


.00006 



Average XFACT = .84719 
Average A(0) = .15414 
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